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oo ■ Abstract 
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0^ . We consider a chemical freeze-out mechanism which is based on a strong medium dependence 

£Nj \ of the rates for inelastic flavor-equilibrating collisions based on the derealization of hadronic wave 



o 



C3 



functions and growing hadronic radii when approaching the chiral restoration. We investigate the 
role of mesonic (pion) and baryonic (nucleon) fluctuations for melting the chiral condensate in 



^ \ the phase diagram in the (T, /i)-plane. We apply the PNJL model beyond mean-field and present 

?H ■ 

an effective generalization of the chiral perturbation theory result which accounts for the medium 
dependence of the pion decay constant while preserving the GMOR relation. We demonstrate 
within a schematic resonance gas model consisting of a variable number of pionic and nucleonic 
degrees of freedom that within the above model a quantitative explanation of the hadonic freeze-out 
curve and its phenomenological conditions can be given. 
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I. INTRODUCTION 



The investigation of the phase diagram of QCD in the plane of temperature T and 
baryochemical potential \Lb is one of the goals of heavy- ion collision experiments, of lattice 
QCD (LQCD) and effective field theory approaches to the nonperturbative sector of QCD. 
Of particular interest are the conditions under which the approximate chiral symmetry 
of the QCD Lagrangian will get restored and whether this transition will be necessarily 
accompanied by the deconfinement of quarks and gluons. More detailed questions to the 
QCD phase diagram concern the order of these phase transitions, their critical exponents 
and fluctuation measures as well as the possible existence of a critical point or even a triple 
point in the phase diagram. Most promising tools for the experimental determination of 
these characteristics are the energy scan programs at CERN, RHIC and at the upcoming 
dedicated facilities of the third generation: FAIR and NICA. The systematic analysis of 
higher moments of distributions of produced particles in their dependence on the collision 
energy and the size of collising systems shall provide answers to the above questions and 
allow direct comparison with predictions from the underlying theory, as provided by LQCD, 
see, e.g., Ref. []J and works cited therein. 

As long as the applicability of LQCD methods is bound to the region of finite T and 
Hb/T <C 1, any predictions for the phase structure of QCD at high baryon densities including 
the possible existence of critical points will rely on effective models. To be relevant for the 
discussion of the above problems these models have to share with QCD the property of 
chiral symmetry and its dynamical breaking as well as a mechanism for confinement and 
deconfinement. 

At present stage, of particular relevance for the discussion of the QCD phase diagram are 
the chemical freeze-out parameters (T? , fi B ) which have been obtained from the statistical 
model analysis of particle yields obtained in heavy ion collisions 

aa. 

One of the most 

striking observations is the systematic behaviour of these parameters with collision energy 
y/s J^-f]], which has recently been given a simple parametric form Q]. It has been observed 
that the resulting freeze-out curve in the phase diagram is closely correlated to the thermo- 
dynamical quantities of the hadron resonance gas described by the statistical model. These 
phenomenological freeze-out conditions make statements about the mean energy per hadron 
(E)/{N) ~ 1.0 GeV, the dimensionless entropy density s/T 3 ~ 7 and a total baryon and 
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antibaryon density % + n§ — 0.12 fm -3 . The freeze-out line provides a lower bound for the 
chiral restoration and deconfinement transition in the phase diagram. Being coincident at 
low densities as inferred from LQCD [8], both transitions need not to occur simultaneously 
at high density thus allowing for an island of a quarkyonic phase [9] between hadron gas 
and quark-gluon plasma with a (pseudo-) triple point [10). 

The question appears for the physical mechanism which governs the chemical freeze-out 
and which determines quantitatively the freeze-out parameters. One aspect is provided by 
the requirement that hadrons should overlap in order to facilitate flavor exchange reactions 
which establish chemical equilibrium. Th is g eometrical picture of freeze-out is sucessfully 
realized in a percolation theory approach Another aspect is the dynamical one: when 
the equation of state (EoS) possess softest points (e.g., due to the dissociation of hadrons 
into their quark and gluon constituents with a sufficient release of binding energy involved) 
which is correlated with the freeze-out curve in the phase diagram, then it is obvious that 
the hadron abundances are characterized by the corresponding T and fiB values 12) . The 
dynamical system got quasi trapped at the softest points for sufficient time to achieve chem- 
ical equilibration before evaporating as a gas of hadron resonances freely streaming to the 
particle detectors. Also the kinetic aspect of fast chemical equilibration was discussed in the 
hadronic gas when accounting for a multi-hadron dynamics 13]. 

All these mechanisms are appealing since they provide an intuitively clear picture, but 
they are flawed by the fact that their relation to fundamental aspects of the QCD phase 
transition like the chiral condensate as an order parameter are not an element of the de- 
scription. 

In the present contribution we develop an approach which relates the geometrical as well 
as the hydrodynamic and kinetic aspects of chemical freeze-out to the medium dependence of 
the chiral condensate. We demonstrate within a beyond-meanfield extension 14H20| of the 
Polyakov NJL model 2lH27| how the excitation of hadronic resonances initiates the melting 
of the chiral condensate which entails a Mott-Anderson type derealization of the hadron 
wave functions and a sudden drop in the relaxation time for flavor equilibration. Already for 
a schematic resonance gas consiting of pions and nucleons with artificially enhaced numbers 
of degrees of freedom we can demonstrate that a kinetic freeze-out condition for the above 
model provides quantitative agreement with the phenomenological freeze-out curve. 
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II. CHEMICAL FREEZE-OUT AND CHIRAL CONDENSATE 



We propose to relate the chemical freeze-out to the chiral condensate in the following 
effective way. As a freeze-out condition for flavor equilibrating reaction kinetics in the 
temperature-chemical potential plane (T, /x) we employ 

r exp (T,/i) = r coll (T,/x) , (1) 

where T exp (T, /x) is the expansion time scale of the hadronic fireball and the inverse of the 
relaxation time for reactive collisions is 

t^i{T,v) = Y^°H n 3 > ( 2 ) 

with i,j = ir, N, ... running over all species in the hadron resonance gas. For the cross 



sections we adopt the geometrical Povh-Hiifner law [28|, |29] 



A(r 2 >(r 2 > , (3) 



where A is a constant of the order of the string tension A ~ 1 GeV/fm = 5 fm 2 . Note that 
this behaviour has been obtained for the quark exchange contribution to hadron-hadron 



cross sections [30 ]. 



A key point of our approach is that the radii of hadrons shall depend on T and /x and shall 
diverge when hadron dissociation (Mott effect) sets in, driven basically by the restoration of 
chiral symmetry. This has quantitatively been studied for the pion [3JJ, where it has been 
shown that close to the Mott transition the chiral perturbation theory corrections can be 
safely neglected and the pion radius is well approximated by 

^(T,/x) = A/- 2 (T, / x). (4) 
It has been demonstrated that the GMOR relation holds out to the chiral phase transition 



where pions would merge the continuum of unbound quark matter [32l. |33| . Since the current 
quark mass is T— and /x— independent and the pion mass is "chirally protected", the T— , 
/x— dependence of the chiral condensate has to be reflected in a similar behaviour of the pion 
decay constant 

fl{T^) = -m (qq) T jMl. (5) 
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The resulting relationship between pion radius and chiral condensate in the medium reads 

The derealization of the pion wave function due to the melting of the chiral condensate as 
expressed in this formula is the most important element of the hadronic freeze-out mechanism 
suggested in this work. 

For the nucleon, we shall assume the radius to consist of two components, a medium 
independent hard core radius r and a pion cloud contribution 

= rl + rl(T^) , (7) 

where from the vacuum values r n = 0.59 fm and = 0.74 fm one gets tq = 0.45 fm. 

For the expansion time scale we adopt a relationship which follows from entropy conser- 
vation, S = s(T,fi) V(r cxp ) = const, and a fireball expansion law V(r exp ). Assuming that 
V(r cxp ) oc r e 3 xp one obtains 

r exp (T, f i) = as- 1/3 (T, t i) , (8) 

with a being a constant of the order one. 

As a first step, we will restrict the discussion to a medium consisting of pions and nucleons 
only, whereby we apply the above relationships. The generalization to a hadron resonance 
gas along these lines is rather straightforward. 

In the following we introduce properties of the chiral condensate in a pion-nucleon and 
in a hadron resonance gas and then applying the freeze-out condition (J2J) we compare our 
model results with phenomenological findings on hadronic freeze-out and its conditions. 



III. CHIRAL CONDENSATE BEYOND MEANFIELD 

We start from the general definition of the chiral condensate 

<*> = as^") ■ < 9 > 

The thermodynamical potential Q(T,fi) can be decomposed into contributions from the 
quark and gluon meanfield and hadronic (quantized) fluctuations in the meson and baryon 
channels 

tt(T, fj) = tt MF (T, fjt) + fi mcson (T, H) + fi b aryon(T, fi) ~ fl(0, 0) . (10) 
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The subtraction of Q(0, 0) removes vacuum divergencies and guarantees that the thermody- 
namical potential, i.e., pressure and energy density, of the vacuum vanish. 
The mean field contribution is given by [36 ] 

fi MF (T,/i) = -2N f J -0-he p + T In [l + 3$e~^^-^ + 3$e" 2/9 ^-^ + e^ 6 *-^] 

}2 
+ + W($;T) , (11) 
4Gr 

where the quasiparticle energy is e p = ^/p 2 + m 2 (T, //), and m(T,fx) = m + a(T,fi). In 
the following we will shorten the notation by dropping the arguments of m(T, //) = m and 
cx(T, //) = a. For the Polyakov-loop potential we take the logarithmic form motivated by 
the SU(3) Haar measure, 



W($,T) 



- a{T) $ 2 + b(T) ln(l - 6 $ 2 + 8 $ 3 - 3 $ 4 ) T 



(12) 



with the corresponding definitions of a(T) and b(T) 16|]. In Ref. fl^ |. a dependence of the To 
parameter on the number of active flavors and on the chemical potential has been suggested. 
Here we will use T = 200 MeV. The mean field value of the traced Polyakov loop is then 
given by <5 = $* = [1 + 2 cos(0 3 /T)]/3. 

The quark mass gap is obtained from the extremum condition dQ/da = 0, equivalent to 



a 

2G 
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7T 



/ dpp 2 ™ [1 - ft ~ ft] 



where the PNJL quark distribution functions are given by 



fi 



(13) 



(14) 



1 + 3$e~^( £ p~^) + 3$e~ 2/3 ( e p~^ + e~ 3/3 ( £ p~^ 
The values of the Polyakov loop $ in the phase diagram are found from a similar gap 
equation corresponding to the solution of the extremum condition dfl/d§ = 0. 
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351 ] for the NJL model 



The contribution of mesonic fluctuations has been discussed in 
and in 0, l^] for the PNJL model. Also recently mesonic fluctuations were included in 
these models within the functional renormalization group approach 
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39]. 



In the following we consider only the on-mass-shell meson contributions and neglect the 
continuum correlations beyond the Mott effect. This is justified because we are actually 
interested in the hadronic freeze-out occurring in a region of the phase diagram which does 
not exceed the limits of the hadronic phase. This contribution is given as 

d 3 k \E M {k) 

M= 



a, 



(15) 



where the index M denotes the actual meson with degeneracy factor du- We restrict this 
sum to the lowest meson state, the pion, with E v (k) = y/W+M*. 

The contribution of pionic fluctuations to the chiral condensate we obtain from 

£ - - ~ AJ £t + mT • (16) 

where the GMOR relation has been used to evaluate dE n (k)/dm Q = M 2 /(2m E 7T (k)) and 
the scalar pion density has been defined as 



=2^J dpP EM ePW - 1 ' 
The last term corresponds to a vacuum contribution which gets eliminated by the vacuum 
subtraction rule ( Imp inherent in the definition of the full condensate (jUj). 

The contribution of baryons to the partition function is considered as an ideal Fermi gas 

fWn(T, li) = - £ dB Jj0f{ ^ + Tln[l + e -«^C*)-«.)] | + ^ ^ -M18) 
For the discussion of the quark condensate we will focus here on the nucleonic contribution 

= = a + , (19) 

where we have used the Feynman-Hellman theorem and introduced the pion-nucleon sigma 
term, ctn = m (dmN / dm ) = 45 MeV [40j. The nucleon scalar density is given by 

n s AT^) = ^- 2 r dpp 2 -^- {f N {T^) + f N {T,-^} , (20) 

where /jv(T, //) = {1 + exp[(A/p 2 + m 2 ^ — yU B )/T]} _1 is the nucleon Fermi distribution, 
tun = 939 MeV is the nucleon mass and /ie = 3/i the baryon chemical potential. The 
vacuum contribution stemming from the zero-point energy term gets removed by the vacuum 
subtraction procedure ( fTUi) . 

From the above discussion, the chiral condensate in a tt — N gas can be obtained as 

/_ v cr M^n S7V {T) a N n sN (T,[i) 
- {qq) = 2G-^^ m • (21) 

Note, that using the chiral limit expression for the scalar pion density n S)7r = d n M n T 2 /12, 
and the GMOR relation (jSJ, we may give this result the form 

T 2 a N n SjN (T,fx) 



(qq) = ( qq ) 



MF 



1 - 



(22) 
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which is known from the chiral perturbation theory 4l|, |42j, but now with a medium- 
dependent pion decay constant. Since the modification of the meanfield contribution stem- 
ming from the quark excitations is proportional to /$, it is now effectively suppressed by 
the Polyakov-loop, compared to the standard NJL model case. 
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FIG. 1: Behavior of the chiral condensate (qq) in the PNJL model at meanfield level (a), with 
pion fluctuations (b), with pion and nucleon fluctuations (c) and for the schematic resonance gas 
medium with 8 pionic and 20 nucleonic degrees of freedom (d). Also shown are the corresponding 
reductions of the chiral condensate relative to its vacuum value for each set of parameters. 



Numerical results for the condensate in the T — \x plane are given in Fig. 1, where we 
display contours of equal condensate and indicate its percentage reduction relative to the 



vacuum value. We use standard values of the NJL model parametrization, e.g., from 
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43|, 



with GA 2 = 2.31825, A = 602.472 MeV, m = 5.27697 MeV, m s>0 = 150 MeV, M w = 140 
MeV, M K = 495 MeV, / w = 92.4 MeV and f K = 93.6 MeV. 

IV. THE FREEZE-OUT CURVE 

In the following we compare predictions of our model on the freeze-out conditions with 
that obtained within the statistical model analysis of particle yields in heavy ion collisions. 

Fig. [2] shows results for different freeze-out curves obtained from the kinetic condition, 
Eq. (jTj) applied to the schematic resonance gas model with variable numbers of pion (d n ) 
and nucleon ((In) degrees of freedom. The results are compared with the phenomenological 
values. From this figure, it is clear that for d n = 8 and djy = 14 ... 20 there is an excellent 



250 




FIG. 2: Freeze-out curves according to the kinetic freeze-out condition, Eq. ([I]), for the schematic 
resonance gas model with different number of pion (d n ) and nucleon (djv) degrees of freedom 
compared to phenomenological values (symbols) [10] . The lines correspond to different choices for 
d w and d/v given in the legend. 

correspondence of our freeze-out criterion with phenomenological values. Increasing re- 
sults in decreasing freeze-out temperatures in the meson-dominated region at small Hb/T. 
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On the other hand, increasing shifts the entire freeze-out curve towards lower T and [1b- 
The freeze-out curves calculated within our model correspond to a reduction of the chiral 
condensate up to 70% of its vacuum value, c.f. Fig. 1(d), so that the hadron-hadron cross 
section in (J2J) according to ([3]) is roughly twice the vacuum value. 

To verify further our model, we evaluate the entropy density s(T,fi) = —dQ/dT, which 
according to the statistical model is related to the freeze-out curve by the phenomenological 
condition, s{T f , fi f B )/T 3 ~ 7. Fig. [3] (left) shows the comparison of the phenomenological 
freeze-out points with the lines of constant s(T,u)/T 3 calculated within our model. The 
dashed line in this figure corresponds to s(T, n)/T 3 = 7, which is the freeze-out line obtained 
within the hadron resonance gas model. 

Fig. E] (right) shows a similar comparison as the left-hand figure but assuming freeze-out 
conditions of constant density of baryons, ub(T^ , + n B (T^, u B ) = const.. The nucleon 
density is obtained from 

nN{T ^ ) = ^l ^ l + e^W-^) ' (23) 

where = 20 is the number of nucleonic degrees of freedom in our schematic resonance gas 
and the antinucleon density, n^(T,fi) = ri]\r(T, —fi). The result of the calculation with the 
above parameters is shown in the right panel of Fig. [3j From this figure one can conclude 
that our effective model reproduces the phenomenological freeze-out line also when imposing 
conditions of fixed total density of baryons. 

V. CONCLUSION AND OUTLOOK 

We have developed chemical freeze-out mechanism which is based on a strong medium 
dependence of the rates for inelastic flavor-equilibrating collisions based on the derealiza- 
tion of hadronic wave functions and growing hadronic radii when approaching the chiral 
restoration. This approach relates the geometrical (percolation) as well as the hydrody- 
namic (softest point) and kinetic (quark exchange) aspects of chemical freeze-out to the 
medium dependence of the chiral condensate. For our model calculations we have employed 
a beyond-mean-field effective extension of the Polyakov NJL model. We could demonstrate 
how the excitation of hadronic resonances initiates the melting of the chiral condensate which 
entails a Mott- Anderson type derealization of the hadron wavefunctions and a sudden drop 
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FIG. 3: Entropy s divided by T 3 (a) and total baryon density n + n (b) as functions of temperature 
T and the baryon chemical potential fiB for the schematic resonance gas model with 8 pion and 20 
nucleon degrees of freedom. The dashed curves represent the phenomenological lines of constant: 
s(Tf,Li f B )/T 3 = 7 in (a) and n B {T f , p f B ) + n B (T f , (J, f B ) = 0.12 fm" 3 in (b). 

in the relaxation time for flavor equilibration. Already for a schematic resonance gas consist- 
ing of pions and nucleons with artificially enhanced numbers of degrees of freedom we could 
demonstrate that the kinetic freeze-out condition can provide quantitative agreement with 
the phenomenological freeze-out curve. The further development of this approach would 
clearly require the replacement of the schematic resonance gas model with the one based on 
the full spectrum of hadronic resonances and detailed studies of their influence on the chiral 
condensate. To this end generalized models for hadron resonance masses are required, 
which reveal their dependence on the current quark mass like, e.g., the model of Leupold 
461 ] . However, it is very likely that with such extension the model results quantifying 
freeze-out conditions will still be consistent with phenomenological findings. 
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